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1 Introduction 

Within last decade there has been achieved large progress in theory of quantum 
computing. Unfortunately the experimental realization of practically valuable 
quantum computers has not been performed mainly due to lack of scalable co- 
herent controllable two-level quantum systems. Semiconductor nanostructures 
were considered as a candidates to basic elements of quantum computer — quan- 
tum bits (qubits). In 1995 A. Barenco, D. Deutsch, A. Ekert, and R. Jozsa 0] for 
the first time offered to use as base states ("0" and "1") of qubits two bottom 
levels of spatial quantization of single-electron quantum dots. To implement 
two-qubit operations it was suggested to use electrical dipole interaction. Soon 
the same authors (A. Ekert and R. Jozsa) Q evaluated coherence of such sys- 
tems and proved that ordinary quantum dots are too incoherent to be quantum 
bits. The variants of such structures were investigated by other authors in pa- 
pers Q . Potential confining at least one of these quantum dots was assumed 
asymmetrical. Distance between the bottom levels was about 10-100 meV. 

Our proposal is to use mesoscopic structures with small (below 1 meV) sep- 
aration between energies of two bottom states. The work frequency of such 
structures, which is proportional to energy separation, is surely to be decreased 
compared to early proposed ones. But processes of spontaneous emission of pho- 
tons and phonons are proportional to polynom (cubic or higher degree) of energy 
of transition. Therefore errors rate per one implemented quantum operation is 
to be smaller compared to common structures. 



2 Structure and principles of operation 

In the proposed structure we offer to use as a qubit a quantum dot with a 
symmetrical potential profile, as shown in Figs. 1, 2. The presence of two minima 
of potential separated by a thick barrier is essential. In a quantum dot there is 
one electron. The presence of the second electron is energetically unprofitable 
because of Coulomb interelectron repulsion. Although the qubit proposed can 
be made in principle of any semiconductor for technological reasons we focus our 
attention on GaAs/AlGaAs structures because this material proved to reliable 
basis for observation of various coherent quantum effects. In GaAs quantum 
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dots at distance between minima r — 10 nm (see Fig. 3) the Coulomb energy is 
about e^/fcr = llmeV, that allows to exclude a spontaneous charging of a dot 
by the second electron. 

To find work frequencies and other working parameters of an offered qubit, 
two-dimensional Schrodinger equation for an electron in a quantum dot in GaAs 
with model potential V (shown on Fig. 2) was solved numerically: 



V = — + Vb cxp 



{wiy 



(1) 



where m = 0.065TOe; I = 20 nm; Vb = 1.5 • lO^^-* J; w = 0.08 0.34 (r = 
11^34nm). 

The obtained wavefunctions of four bottom states of an electron for this 
potential are shown in Figs. 4-7. Two bottom energy levels of an electron in 
such point correspond to wavefunctions ^'i (Fig. 4) and ^'2 (Fig. 5). For logic 
"0" and "1" it is convenient to take the normalized sum and difference of 4'i 
and \E'2 rater than states with the certain energy: 



|0) = (Vl/i + VE'2) , |l) = i=(vI/i-*2). 



(2) 



The given states correspond to almost complete localization of an electron in 
one of minima of potential, as shown in a Figs. 8, 9. It allows to implement write 
of the input data and reading of results by methods of single-electronics with 
the help of the read-out gates, as shown in a Fig. 1. The central gate (control 
gate) serves for downturn of a potential barrier while implementing quantum 
unitary transformations, as will be shown below. 

2.1 Implementation of one-qubit unitary operations 

In the beginning qubit is not affected by operations and it is in the state 



*(0)==co|0)+ci|l). 
Electron wavefunction evolve according to Schrodinger equation 



where Hamiltonian 



2m 



(3) 



(4) 



(5) 



where V is given by 

Different voltages on the control gate correspond to different barrier heights 
Vb- Therefore energies of the ground and first excited state also depend on 
Vb- (?iwi,2 = ^^^^1.2 (Vb)) ■ When barrier is high, the energy levels of two 
bottom states practically merge so two bottom levels evolve with the common 
frequency w: 
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vl/(t) = e— *(co|0)+ci|l)). (6) 

The downturn of a barrier results in an inequality of frequencies (Aw = 
uj2 — ^i > 0) and to periodic rotation of vector of state of qubit in basis (|0) , |1)). 

^(t) - ^ (CO + ci) vPie— 1* + 71 - ^26-^'^^*, (7) 

that is 



^r(i) = e''^'^'+'^'-^'^^[{cQCOsAujt/2 + icismAujt/2)\0) + (8) 
+ (ci cosAujt/2 + ico sinAwi/2) |1)]. 

Having applied a pulse of a positive voltage of the certain duration tjsfOT, 
equal to tt/Aw, we (an insignificant common phase factor exp [—inuji/ {lo2 — wi)] 
can be neglected) shall transfer a state of a qubit (^) into a state (as shown in 
Fig. 10): 

NOT(*o) = ci|0)+co|l). (9) 

So, with the help of the given procedure it is possible to exchange amplitudes 
at "0" and "1", that is to carry out unitary NOT operation. Changing the 
duration of a pulse it is possible to carry out rotation of qubit state to any 
required angle. 



2.2 CNOT gate implementation 

For construction of the universal quantum computer it is also necessary to be 
able to realize at least one nontrivial (not decomposable into a sequence of one- 
qubit gates and permutations) two-qubit operation. We consider realization of 
CNOT operation. (Depending on state of control qubit, on target qubit should 
be carried out either operation of identity or NOT) . For implementation of the 
CNOT gate between the neighbour qubits the Coulomb interaction is used. 
We arrange two qubits as shown in a Fig. 3. To exclude exchange effects it 
is supposed that the qubits are separated by completely opaque barrier. Then 
Hamiltonian of electron in target qubit becomes 

Ht = - — A + V + Vc, (10) 
Im 

where Coulomb potential due to electron in control qubit 

Vc{.,y)= lidudv ^ \''cM\'e- 

K^{x^vf + {y + R-uf 

where {u^v) is a wavefunction of electron in control qubit, u,v are coordi- 
nates in coordinate system of control qubit, e is the electron charge and i? is a 
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separation between qubits' centers. Height of barrier dividing the right (target) 
qubit depends on state of the left (control) qubit. The addition to the height 
of barrier of the target qubit can be roughly estimated in the following way: 
when control qubit is in state |0) or |1) we can assume that electron is a point 
charge located in electron density maximum. Then the Coulomb potential in 
the barrier dividing target qubit will be 

Vc{0,v)^e^/n{y + R + sr/'i), (12) 

where s corresponds to state of control qubit: s = 1, when control qubit is in 
state |0), and s = — 1, when control qubit is in state |1). 

In such structure CNOT gate can be implemented as follows. Having slightly 
opened the barrier in a target qubit with the help of its control gate it is possible 
to achieve application to the target qubit the operation of identity, when the 
control qubit is in a state |0) and operation NOT, when the control qubit is in 
a state |1), that is operation CNOT. Consider a case, when the control qubit is 
in one of base states. Because of Coulomb influence of a control qubit on height 
of a barrier in a target qubit the duration of the NOT operation in these cases 
will differ {Inoto i jvoti ) • Duration of a pulse on a gate of a target qubit we 
choose in the following way: 

tcNOT = %OT0 tNOTl/{tNOTl — %OTo)- (13) 

Varying amplitude of a pulse we shall achieve, that the ratio 

tNOTl/^{tNOTl — tNOTo) (14) 

will be an integer (it is the large value, as in our case interaction is weak, so by 
small change of amplitude we attain the nearest integer value) 

tjVOTl/2(iArOTl — tNOTo) = N. (15) 

If the control qubit is in state " 0" , then the action of a pulse is equivalent to 
consecutive application of even {2N) of number of operations NOT to the second 
qubit, that is operation of identity. If the control qubit is in state '' 1" , then the 
action of a pulse is equivalent to consecutive application of odd {2N — 1) number 
of operations NOT to the second qubit, that is NOT operation. 

While modelling the interaction was calculated directly by the Coulomb's 
law (in all cells of a two-dimensional grid the field created by a partial charge 
of a control qubit form all cells) was calculated. Exchange effects were ne- 
glected. The dependencies of minimal durations of operations NOT and CNOT 
depending on geometrical parameters are shown in Fig. 11, 12. 

3 Numerical modeling 

For calculation the method of simple iterations was used. Potential should be 
symmetrical and have 2 minima, that means presence of a barrier between these 
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minima. Similar potential for a two-electronic quantum dot was used in [gj : he 
has offered to use polynomial potential of the 4-th order like 

Vi.,y)^^(^^{x'~a'f + r/'^ (16) 

with minima in points ±a. 

In this work as potential of a quantum dot with the built-in tunnel barrier 
was used potential (|l]). Such potential can be varied by changing the following 
parameters: I (characteristic size of a qubit (~ 20 nm)), Vb (height of a barrier 
(~ 1 eV)) and w (the width of a barrier in relation to I (~ 0.1)). Varying values 
of these parameters, it is possible to tune the characteristics of structure. 

Consider an electron strongly limited in a direction, perpendicular to a sur- 
face of a heterojunction. Thus distances between levels of spatial quantization 
(2D electron gas subbands) in this direction is about 100 meV. At low temper- 
atures and weak influences the electron remains in the bottom subband. Thus 
its motion can be considered as two-dimensional, and the wavefunction of an 
electron is factorized into 

<^{x,y,z)^^{x,y)g{z). (17) 

For wavefunction (x,y) we have solved 2D Schrodinger equation with the 
potential V using simple iterations method with orthogonalization on each iter- 
ation. We have modelled an area 60 x 40 nm^ using mesh with the same step in 
both directions 0.5 nm. We have iterated 4 bottom state wavefunctions control- 
ling their orthogonality. The mean number of iterations was about 10000. To 
find energies of states of an electron, the mean value of an energy for a wavefunc- 
tion 4" under the formula {^\H \ where H is the operator of Hamilton with 
real potential was calculated. After that the program show us the diagram of 
the bottom part of a spectrum of system, and also frequency of transitions from 
the exited energy levels on basic. These frequencies were also used in research 
of processes of decoherence in structure. Besides half of time of transition from 
the first exited level on basic and back also was accepted for the minimal time 
of operation of the NOT gate. 

For the nontrivial two-qubit CNOT gate it is necessary to include interaction 
between control and target qubits. In studied structure the electrical (mostly 
dipole) interaction is used. For this purpose we arrange a target qubit about 
an end face of a control qubit, having turned it on 90 degrees in a plane of 
structure. At such arrangement a state of an electron in a control qubit (the 
electron is farther or closer to the middle of a target qubit) effectively reduces 
or increases height of a barrier in a target qubit, that changes time of operation 
of the NOT gate. As the points of localization of an electron in a target qubit 
are located symmetrically in relation to the control qubit, so the change of a 
state of a target qubit does not render influence in the first order on a state of 
a control qubit. 

For calculation of interaction it is necessary to calculate potential created 
by an electron of a control qubit in the field of a target qubit. For this purpose 
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it is possible to take advantage of two methods: solution of a Poisson equation 
or calculation of an electrostatic field using Coulomb's law. The second method 
was used. For this purpose it was supposed, that in each node of a mesh in the 
area of a control qubit there is a point charge, which value is equal to value of a 
charge of an electron multiplied on a value of function of spatial distribution of 
an electron at this node of a mesh. As function of distribution it is natural to 
take a square of the module of a wavefunction. As the function of distribution 
is normalized, this operation simply represents an electron by system of point 
charges with a cumulative charge equal to a charge of an electron. 

The durations of the NOT gate for two extreme states of an electron in a 
control qubit were obtained: when the electron wholly is in most distant from 
a target qubit and when the electron wholly is near to a target qubit. Knowing 
these times, it is possible to calculate duration of operation of the CNOT gate. 

4 Coherence of structure 

Significant difficulty interfering creation of the large scale quantum computer is 
the problem of decoherence of a quantum state because of interaction with an 
environment bringing in errors. It is proved that if the decoherence occurs 
slowly enough, that, in particular, means, that at calculation occurs no more 
than i5 failures for one computing step (by various estimations 6 should lie in a 
range from 10~^ to 10~^), with the help of special error correcting algorithms 
and codes demanding polynomial increase of computing expenses, modeling 
functioning of ideal (coherent) quantum computer and steady implementation 
of any quantum algorithm is possible. 

While calculating the decoherence in offered structure the low-temperature 
limit (T — i- 0) was considered. The given approximation is justified, as the mod- 
ern cryogenic engineering allows to carry out functioning nanoelectronic struc- 
tures at temperatures down to several millikelvins, that is sufficiently lower 
than the distance between bottom and first excited levels. The case of high 
temperatures represents only academic interest because of inevitable fastest de- 
coherence and impossibility of correct work of the quantum computer. However, 
in solid-state structures even at absolute zero of temperature the processes of 
decoherence owing to spontaneous emission of photons or acoustic phonons with 
transition of an electron from excited to the basic level are possible. These pro- 
cesses will determine the degree of coherence in our structure. We separately 
investigated spontaneous emission of photon, of deformation acoustic phonon 
and of piezoelectric acoustic phonon. 

4.1 Emission of photon 

Consider process of decoherence through emission of photon. The transition 
from first exited on the basic level is dipole. The probability of dipole transition 
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is given by the known formula [Q 

W..^^M^ (18) 

where the dipole moment 

d = y"d2r$'(r)*er$(r). (19) 

From the symmetry of wavefunctions it follows that of a component dy will be 
zero, and x-component can be calculated through integral not on all space, but 
only on half-space x > 0: 

d^^2 J dxdy^' (r)* ex^{r). (20) 



x>0 



As the integrand at large r exponentially decreases with increasing of x, integral 
can approximately be evaluated, having reduced integration on half-space to 
integration on area < a; < Xmax ■ To estimate this integral we can replace x in 
an integrand with its maximal value, i.e. 

dx < 2ex„iax J dx dy $' (r)* $ (r), (21) 

0<a;<2;max 

and as the wavefunctions on area of integration practically coincide, integral 
will be equal to 1/2, as the integration is made on half-space, and integral on 
all space from a condition of a normalization of wavefunctions is equal to 1. In 
our case of wavefunctions have on a; > maxima in some point r, behind which 
exponentially fall down, and as XmO-x it is enough to choose 2r. Then d^ < 2er, 
hence estimation for probability of spontaneous emission of photon to look like 



4.2 Emission of acoustic phonon 

The probability of relaxation from the excited state to the ground one, by 
emission of acoustic phonon, is calculated by formula: 

wio = Y \^'^\^ ^ - hsq) , (23) 

where M — {f\T \i) is the initial state (the electron is in the excited state 
with an energy £i), |/) is the final state (the electron is in the base state with 
an energy Eq, Eiq = Ei ~ Sq, the phonon carried away an energy hsq), T is the 
transition operator) is the matrix element of transition appropriate to emission 
of phonon, s is the speed of sound (5.2 • 10'' m/s in GaAs), q is a wave vector 
of the phonon. 
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Disturbance created by one phonon sq 

1 



'^^^ - 13/2 



3 1 1/2 

e'i"-t;,q (r) , (24) 



_2MoU!sq_ 

and 

z;,q(r) = ^K„(r)d?qe^'i('^-'-), (25) 

aa 

where Mq is the mass of an elementary ceU, Oq is the volume of an elementary 
cell, c? is a unitless vector of polarization orthonormalized by the condition 

(d?q) * d^,q, = Afo<5,q,,-q. , (26) 

a 

where Ma is a mass of atom a. From a normalization follows, that for acoustic 
phonons at q ^ all atoms in a cell are displaced equally: 

d^q = d,q, |d,J = l. (27) 

It is convenient to define a total probability of emission of phonon, i.e. prob- 
ability to emit though any phonon, considering, that the crystal is limited to 
a normalizing volume , and the allowable values of a pulse are quasidiscrete. 
Then the probability of transition is probability that in unit of time electron will 
make a transition from one quasidiscrete state in another; the dimensional rep- 
resentation of it is . The certain thus probability depends on a normalizing 
volume. Probability of transition of an electron from the given state 

W = Yw^. (28) 
q 

As the allowable values of q are located very richly, on distance 2?i/L, it is 
possible to pass from a summation on q to integration by a rule 

Then 

«. = .= /|2,.„. (30, 

4.3 Emission of piezoelectric acoustic phonon 

In 1961 A.R. Hutson considered interaction of electrons and acoustic waves due 
to piezoelectric effect ^ . The matrix element of interaction of an electron with 
an electromagnetic field is given by the formula 

MpA= I d^r^' [r)* eip{r)^{r), (31) 
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where $ (r) is wavefunction of an initial state of electron $' (r) - of the final 
state \f) , (fi (r) is the macrofield created by a phonon. It can be found from the 
Poisson equation 

VV = 47rdivP, (32) 

where P is a dipole moment of unit of volume arisen at deformation of lattice. 
Use of the Poisson equation instead of complete system of the equations of the 
Maxwell corresponds to the indefinitely large speed of light c. This assumption is 
justified, as c is much greater of phase speed uj/qol phonons, that are interesting 
for us. 

At homogeneous acoustic deformation 

Pj = PjkiUki, (33) 

where /? - piezoelectric constants. Having substituted ( |33| ) into (^), it is possible 
to find a field (p. If we are interested in a field ipsq, created by single phonon 
sq, from the Poisson equation we have 

ips^ = -i(47r/g2)qp^^, (34) 

where P^q is a polarization created by one phonon sq. 
As a result for acoustic phonons it turns out B 



/ ^ 1 



1/2 

e^'i'-""-'Aq + c.c., (35) 



_2Mqujs 

where the effective piezoelectric constant of wave sq is entered: 

/3sq = 47rejefe/3ifcjdiq, e^q/q; (36) 

it depends only on a direction of phonon propagation and from a polarization. 
Thus, the probability to emit a piezoelectric acoustic phonon 

wio = (e/3sq)^ |/ (q)l^ 5 (eio - Tisq) , (37) 

where 

j dV$'(r)*$(r)e^(9"=^+«''?'). (38) 
With the account of for a total probability we finally have 

= TTT^ / (eAq)' 1^ (q)l' -5 (eio - hsq). (39) 

(27r) peio J 

For calculation of probability it is convenient to pass to spherical coordinates. 
This transition was carried out by a rule 

q-x ~ q cos 9 cos ip, 

qy = q cos 9 sin ip, (40) 
qz = q sin p. 
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In spherical coordinates the expression for calculation of probability will be 



WPA^-^ d9d^dqq^cose{el3s^f\I{ci)\^ 5{ew-hsq) (41) 

{2tt) pew J 

We get rid of delta-function under integral, having made integration on dq: 

WPA= f (t^)'f- / d9 dv COS 0{e(3s^f\I{<i)\\ (42) 
and finally we get 

= . /'!2 , / d0d^ coseiePs^f (43) 
Sn'^pn s-' J 

Vector of displacement of atoms is 



(44) 



We consider cubic crystals of symmetry classes Oh and T^. In crystals Oh 
with the center of inversion (for example, in Si, as it is possible to show from 
reasons of symmetry, P = 0, i.e. the homogeneous deformation does not create 
macrofields. In crystals Td without the center of inversion (for example, GaAs) 
the tensor (3ikj has only those components, in which all three indexes i, k,j are 
various, and all these components are equal. Therefore 

47r/3 

/3q = (qxQyC + qyqza + qzqxh) , (45) 

where (3 is constant (ei4/Ko, where ei4 is sole piezoelectric constant of cubic 
crystal (0.16 C/m^ for GaAs [^[l^), kq - dielectric permeability (12.8 for GaAs 
|, 10|)). In next sections we will fulfil calculations of probability of emission of 




piezoelectric acoustic phonon for GaAs. 
4.3.1 Transverse phonons 

The phonon with a transverse polarization propagated in any direction, can 
be decomposed on basis consisting of two polarizations, perpendicular to each 
other. Choose the first direction of polarization not causing displacement of 
atoms in direction z. The vector of displacement of atoms for the first direction 
of polarization also should be perpendicular to the wave vector q and to be 
normalized: 

1 / '^y 

dri = , -fe I • (46) 

+ 
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For such wave an effective piezoelectric constant 

Priq = ^— V cos^ 6 cos 2(p sin ( 



(47) 



= — 477/3 V cos^ 9 cos 2(fi sin 6, 
and under integral there will be an expression 

{e^Tidf = (47re/?)%os2 6lcos^ 2<y!?sin2 61. 



(48) 



The vector of displacement for the second direction of polarization should 
be perpendicular to the wave vector q, to the vector d and to be normalized: 



dT2 = 



9; 



1 \l li+'li 



\ 



(49) 



Effective piezoelectric constant for this direction of polarization 

/3T2q = ^— (cos^ + 3cos3^) cost/jy'sin^^ = 

= —77,9 (cos ^ + 3 cos 30) cos (/j-y/sin^, 
and under integral there will be an expression 

(e/?T2q)^ = (TTf /3) ^ (cos + 3 cos 30)^ cos^ sin^ i-p. 



(50) 



(51) 



Summarizing on polarizations of phonons, thus, in the integrand we obtain 
an expression 



(e/JTc)' = (eW + (e/3T2q)'=(^ 



cos X 



(52) 



X (4 (7 - 9 cos 261) cos 4(^ cos^ 61-28 cos 20 + 9 cos 40 + 27) . 

Finally, the total probability of emission of phonon with a transverse polar- 
ization is given by the following expression: 



Wtpa 



eio {epy 



1 dOdif 


J COS0 \ 









cos y X 



X (4 (7 - 9 cos 20) cos 4ip cos^ 0-28 cos 20 + 9 cos 40 + 27) . 



(53) 
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4.3.2 Longitudinal phonons 



The phonon with a longitudinal polarization propagated in some direction, will 
cause displacement of atoms in the same direction: 



1 



IL = 



qx 

Qy 



(54) 



Effective piezoelectric constant for such vector 

PLq = -^SqxqyQz, 
and the integral expression will get 

\ 2 



(55) 



(e/3Lq) 



(IW)^ , .2 

^ — (qxqyqzj 



[Uirepf cos'^ 9 sin^ cos^ ip sin^ (p. (56) 



Finally, the total probability of emission of phonon with a longitudinal po- 
larization is given by the following expression: 



W, 



LPA 



p{hs) 



do dip 



I £10- 



cos 9 



hs ' 



cos^ 9 sin^ 9 cos^ p sin^ p 



(57) 



4.3.3 Emission of deformation acoustic phonon 

The matrix element of interaction of an electron with a deformation field intro- 
duced in 1950 by J. Bardeen and W. Shockley [0 is given by 

Mda = J d^r^'{r)*w{r)^{r), (58) 

where value w{r) is called deformation potential. The homogeneous deforma- 
tion of crystal due to long-wave acoustic phonons is described by a tensor of 
deformation Uij. Therefore in the lowest order on displacement of atoms defor- 
mation potential of acoustic phonons can be expand into a series of Uij and to 
write down 

w = (59) 

(the summation on repeating indexes is performed). Values S*^ (of the units 
of energy) are referred as constants of deformation potential. As the tensor u.y 
is symmetrical, the same property has the tensor of constants of deformation 
potential S*-' . The number of independent constants is determined by symmetry 
of Brillouin zone in that point k, in which the influence of deformation on a 
spectrum is studied. 
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The elementary case takes place in a cubic crystal, when the extremum of 
the band is located in a point k = 0. The symmetrical tensor of the second rank 
in this case is reduced to one constant 



T^jj _ x. t: 



(60) 



Then 



W — ^M, 



U = Uii + U22 + M33, 



(61) 



where u - relative change of volume as a result of deformation. In this case de- 
formation of shift which is not changing of volume, does not result in occurrence 
of deformation potential. From ( p4[ ) and ( |25| ) one obtain 



1 



1/2 



(62) 



Substituting 
phonon 



into (|5^), one find deformation potential from one acoustic 



1 



L3/2 



1/2 



e^'""-'"-*S,qZg + c.c., 



where the effective constant of deformation potential of wave is included 



1 



"^2^ Kq*- +rfiq*) , 



(63) 



(64) 



it depends only on a direction and polarization of phonon. In our case of cubic 
crystal we have S^q = S (7 eV for GaAs for longitudinal phonons and Ssq = 
for transverse ones |l^ ). 

So, in Si and GaAs the deformation acoustic phonons only with a longitu- 
dinal polarization can be emitted. The matrix element will finally look like 




Mda = «5 



Then probability of deformation acoustic phonon sq emission 



(65) 



WW 



psL- 



r9|/(q)| (5(eio-?isg) 



(66) 



Finally, substituting (|66|) into (|3C|), passing to spherical coordinates by the 
rule ( ^o|) and making integration over delta-function, we have a total probability 
of generation of deformation acoustic phonon 



2^3 



2tt 



"/2 



Wda 



10 



87r2rs5 



dip 



dO cos 9 



cost' 

I I £10^ — ,</5 
ns 



(67) 
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4.4 Results and discussion 

The results of calculation of processes of decoherence by structure owing to 
spontaneous emission of particles are given in a Fig. 13. As it can be seen 
from the graphics the prevailing mechanism of decoherence is the emission of 
polarizing acoustic phonons. However, even this process has small probability 
at durations of step appropriate to a GHz range of clock frequencies for wide 
structures (wider than 15 nm). The energy splitting versus distance between 
electron density maxima r shown in Fig. 14 is far above temperature limit 
(5 mK). It is worth of noting, that the process of emission of polarizing acoustic 
phonons is characteristic for materials such as GaAs, in Si, for example, owing 
to symmetry of lattice they will not be generated. Thus, in Si the prevailing 
mechanism of decoherence will be emission of deformation acoustic phonons. 

5 Conclusions 

New quantum bit is offered on the basis of spatial states of electrons in sym- 
metrical semi-conductor quantum dots controlled with the help of voltage on 
electrodes. The quantum-mechanical calculation of states of the offered qubit is 
fulfilled. The operation of the one-qubit gate which is carrying out unitary trans- 
formation, and not trivial two-qubit CNOT gate are simulated. The processes 
of decoherence were investigated. Our study comprises spontaneous emission of 
photon, deformation acoustic phonon and piezoelectric acoustic phonon. The 
work frequencies of the quantum register being built using offered quantum bits, 
lie in a range convenient for electronic control and achieve 1 GHz. The offered 
quantum register is scalable and its size is not limited. The structure reveals 
sufficient degree of coherence that allows by using of the appropriate methods 
of error correction to work during unlimited time. It is necessary to note, that 
offered quantum bit can be realized physically at an existing level of cryogenic 
engineering and nanoelectronic technology. 
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Figure captions 



Fig. 1. Quantum dot qubit. 

Fig. 2. Quantum dot potential profile (grayscale, whiter denotes the higher 

potential). 

Fig. 3. Sketch of CNOT gate (r is a separation between electron density 
maxima; i? is a separation between qubits' centers) 

Fig. 4 7. Wavefunctions of four bottom states of an electron. 
Fig. 8. Wavefunction of |0) qubit state. 
Fig. 9. Wavefunction of |1) qubit state. 

Fig. 10. NOT gate operation time diagram. is a control gate voltage, Alo 
is an energy separation between two bottom states. 

Fig. 11. NOT gate duration versus distance between electron density max- 
ima r. 

Fig. 12. CNOT gate duration versus distance between dots centers R at 
different r. 

Fig. 13. Decoherence due to spontaneous emission of photons and acoustic 
phonons versus r. 

Fig. 14. Energy splitting sio versus distance between electron density max- 
ima r. 
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